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Abstract 



A simple backreaction problem in quantum mechanics, the full quantum 
anharmonic oscillator, and quantum parametric resonance are studied using 
Renormalization Group techniques for global asymptotic analysis. In this 
short note this technique is adapted for the first time to operator problems. 
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I. INTRODUCTION AND MOTIVATIONS 



In many different contexts perturbation theory fails miserably because of the growth 
of higher order terms, contrary to the basic perturbative assumption. This secularity is 
present in both classical and quantum theories, and pervades the motivation for the search 
for analytical methods to improve on perturbative expansions. 

There is a set of problems of particular interest to us where improvement on perturbative 
expansions has shown itself to be of immediate necessity, and of direct physical meaning, 
namely, (p)reheating in inflationary cosmology We shall not deal directly with this 

problem, though, and we will content ourselves with an analysis of two simple quantum me- 
chanical problems. The choice of problems and treatment thereof, however, will be inspired 
by the mentioned physical setting. 

The method we use is novel, and subsumes many other previously known ones. It is the 
Renormalization Group (RG) method for global asymptotic analysis, as advocated by Gold- 
enfeld and collaborators [§,§J, which we extend to operator problems in this short note. The 
key idea of the RG method for global asymptotic analysis is the introduction of a time pa- 
rameter, additional to the initial value point, in such a way that the perturbation expansion 
is valid in the vicinity of the introduced time parameter. The coupling constants/constants 
of motion/initial conditions (depending on your viewpoint and background one or another 
of these descriptions will be more suitable) are turned into running constants, that is to say 
that these constants are suitably modified by the change of the introduced time parameter. 
On the other hand, the solution itself cannot depend on the additional, new time parameter, 
so derivation with respect to the latter of the perturbative solution will impose evolution 
equations for the running constants. These equations are then solved for the running con- 
stants, and on substitution in the perturbative expansion, together with the choice that the 
time parameter is constantly updated to be time itself, we obtain an improved solution. 

This method has the clear advantage over multiple scales perturbation analysis that 
no a priori determination of the scales that appear in the problem is necessary, and a naive 
perturbation expansion is enough as a starting point. Many examples and illustrations of 
this advantage of the RG method can be found in the works of Goldenfeld and collaborators. 

The multiple scales method itself has been applied to operator problems PJlOfl, but 



not the RG method, and this paper is the first example, to the best of our knowledge, of 
the application of the RG method to operator problems. We shall choose the quantum 
anharmonic oscillator and the phenomenon of quantum parametric resonance as our case 
studies, because of the important role they have traditionally had as theoretical laboratories 
for new perturbative methods, and because of their paradigmatic character in the context 
of cosmological (p)reheating. Our results are comparable to all other methods, and since 
they are obtained with a modest effort, we think that the RG method is highly competitive 
in the operator context as well. There are a number of areas where its usefulness might be 
proved, such as quantum optics, but that we leave for further work. 

II. SETTING OF THE PROBLEM 

The problem we shall first address is the back-reaction problem of the quantum an- 
harmonic oscillator; that is to say, the problem posed by the backreaction of quantum 
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fluctuations on the expected value of the position. Consider thus a lagrangian of the form 



-mq 2 — -mu) 2 q 2 — -Amg 4 . (2-1) 



The equation of motion in Heisenberg's picture reads 

d 2 q 



df2 +uj 2 q + Xq 3 = J (2.2) 

which is, of course, an operator equation, the quantum Duffing equation. Let us now de- 
compose the q operator into a c-number times the unit operator and a fluctuation operator, 
in such a way that the expectation value of q is equal to the mentioned c-number: 

q = if + i, (2.3) 

where <p = (q), (£) = 0, the unit operator is omitted, and (•) stands for the expectation 
value of an operator over a given state, which, as we are in Heisenberg's picture, is a time 
independent state, and carries the initial value data for the operator evolution equation. 
Both the c-number ip and the operator £ are time dependent, of course. 
Taking the expectation value of the operator equation (|2.2|) we obtain 

^ + ^V + A(^ 3 + 3^ 2 > + <a)=0, (2.4) 
and using this result, together with ( |2.2j ), we can write 



y + (u 2 + 3 v) e = sa^ ((e) - e) + a (<a - e) . (2.5) 



It is clear that the pair of equations ( [2.4D and (|2.5| ) is equivalent to the operator equation 



(2.2) together with an initial condition provided by the state over which the expectation 



values are computed. 

So far we have followed the computations in Ehrenfest's theorem [11], and the result is, 
as is well known in elementary quantum mechanics, that the evolution of the expectation 
value of the q operator is not given by the classical evolution equations: there are quantum 
corrections because the potential is not quadratic. 

The customary way of handling this kind of problem is to revert to the interaction picture 
and use time independent perturbation theory. We shall depart from this well-trodden route 
and study the large time asymptotics of the expectation value from a truncation of the exact 
equations ( |2.4| ) and (|2.5| ). The truncation, which we will call the one-loop approximation 
because of the analogy to the one-loop truncation in quantum field theory, leaves us the 
following coupled equations, one c-number and the other operator valued: 



dV 
dt 2 
d 

dt 2 



+ (^ 2 + 3A(a)y = -A^ 3 , 
^| + (u; 2 + 3A</)£ = 0. (2.6) 



3 



Notice that for these truncated equations the equation for fluctuations is linear, with a time- 
dependent frequency related to tp 2 . We could now try a naive perturbation expansion for 
these equations, of the form 

ip = ip + \<px + 0(A 2 ) , 
£ = £o + A£i + 0(A 2 ), 



which, when inserted in eqns (|2.6|) , yields 

d 2 V?o 2 
— — + w ¥?o = 



dt 2 

1 i , ,2,„ Q/t2\,„ ,„3 



dV 
dt 2 

d 2 6 



dt 2 



d 2 ^ 



! 1 w 2 fi = -3#&,, (2.7) 



dt 2 



etc. It now becomes apparent that even within this simple truncation the perturbative 
expansion cannot work, since secular terms will appear and the perturbative hypothesis, 
namely, that the terms in higher order in A are smaller than the previous ones, will stop 
holding after a time of the order 1/A. This is the issue we shall now address by means of 
the RG method. 



III. APPLICATION OF THE RG METHOD 

The zeroth order equations in the series ( |2.7| ) are solved by 

ipo = Rcos(ut + 6) , 

£ = I (a ] e wt + ae- wt ) , (3.1) 

where a and are formally adjoint to each other. The commutation relations among these 
two operators are fixed by imposing canonical commutation relations for q and its conjugate 
operator p. We obtain 



and, therefore, choosing I = yh/(2mu), a and can be understood as annihilation and 
creation operators for a harmonic oscillator, implying that we have a way of defining a 
number operator iV = o^a that will provide us with information about the magnitude of 
the fluctuations. 

So far everything has been straightforward. Problems first arise in computing the first 
order corrections, since secular terms will be generated. Consider a time r for which the 
first order corrections vanish. The solutions to the first order equations will be 
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<Pi = ^(t - r) x [e^ (i?V* + Al 2 ((a^)e^ + 8 / 2 ((TV) + I) e «) - 

- e-™ (R 2 e- ie + Al 2 (a 2 )e id + 8l 2 ((N) + ^) e"*) J + r.t. , (3.3) 



and 



& = _ r ) ( e -* (2«t + ae 2ie ) - e"*"* (2a + aV^)) + r.t. , (3.4) 

where r.t. stands for regular terms, that is to say those which are bounded when t — r tends 
to infinity. 

It is clear that for small t — r the perturbative expansion has no a priori reason not 
to hold. Now, the crucial point is to realize that the choice of r corresponds to a choice 
of initial time for perturbations, which, in other words, means that a particular value for 
the magnitudes R and 9 and a particular choice for the (constant) operators a and a^ has 
been chosen. If we were to allow those to change with r, we would be able to readjust the 
perturbative series as we went away from a given r for a fixed t to a different r', around 
which the perturbative expansion for t — r' would now proceed. The change in R, 9, a and 
cr would compensate the change in r. 

On the other hand, the solution to the differential equations which we intend to solve 
by a perturbative expansion cannot depend on r. Once initial values R and 9 and initial 
operators a and a^ are given, the evolution of (p and £ is fixed. Therefore, we obtain the 
first order RG equations by allowing the stated dependence of the (number and operator) 
parameters in r, deriving the expressions ipo + Xipi and £o + <^£i with respect to r, and setting 
those derivatives equal to zero. To keep consistency with the perturbative expansion, the 
terms of the form XdR/dr and similar are discarded, since they are of order A 2 , which is 
being disregarded. 

Reordering of the first order RG equations and separation of the coefficients of e tuJt and 
e~ tuJt , together with some manipulation of the operator equations, lead us to the following 
closed set of five coupled differential equations for five ordinary functions of r: 



1 dR 
Rd^ ~ 


3iA/ 2 / 

4c y {[ 


d9 

d7 " 




d(a 2 ) 


-3i\R 2 


dr 


8lu 


d((«t) 2 ) 


3iXR 2 ( 


dr 


8uj V 


d(N) 


SiXR 2 


dr 


8u 



(4(a 2 ) + 2e- 2W (N) + e-™) 



(JY)e-™ - (a 2 )e 2ie ) (3.5) 

These equations, after a little reorganising in order to write them as a set of real equations 
with no free parameters, can be solved numerically without much further ado, leading for 
predictions for the long time behaviour of R(r) and 9{t). 
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Consider now that we constantly update r to be t. The secular terms will disappear in 
if i, and their effect will show up through R and 9, now understood as functions of t: 



if = R(t) cos(ut + 6(t)) + e(r.t.) , (3.6) 

where r.t. stands for regular terms, i.e., well controlled small perturbations to the first term. 
Moreover, it can be seen directly, from the analysis of equations (|3.5|), that 



dR l2 d(N) . . 

R— + 2l 2 ^-^ = 0, 3.7 
dr dr 

which is the equation that implements the conservation of energy to this order. 

The numerical solution for the initial conditions y?(0) = ^8u/(3X), (N) = 0, 6(0) = 0, 

(a 2 ) = and ((at) ) = is as shown in figure [I], where (pit) is depicted. The units used are 
arbitrary. 



IV. DIRECT APPLICATION OF THE RG METHOD TO THE QUANTUM 

DUFFING EQUATION 

In the previous section we have demonstrated the usefulness of the RG method in elim- 
inating the secular terms for the truncated equations obtained at one loop. We shall now 
consider the whole quantum Duffing equation, eqn. ( |2.2|) , and apply the same techniques 
for the complete problem. 

Let us perform a simple perturbation expansion in the A coupling constant, q = qo + 
Xqi + 0(X 2 ). The solution to order is simply 



qo 



h 'pe^ + fte iuJt ) . (4.1) 



2mw 



We have written {3 and (3 in order to make the distinction with the operators cr and a, which 
were creation and annihilation operators for fluctuations in the previous section, whereas ft 
and f3 are creation and annihilation operators for the full oscillator problem. The first order 
equation presents resonance and, therefore, secular terms. Writing down just the singular 
(i.e., secular) part of qi, which we call gi jS , 



+ 0$ + ftp 2 ) e- iu}t ) . (4.2) 



Let us now allow a dependence in r of the pair of operators j3 and ft . Imposing the RG 
condition dg/dr = 0, we obtain 

% + ^ {ftP' + + Pft) = 0{X 2 ) , (4.3) 

and the hermitian conjugate thereof. We notice that Af = ft (3 and [f3, ft] are constants 
under the flow of r, which allows us to solve these equations in the form 
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/3t( T ) = e 3iAWVTr/(W)^t(o) , (4.4) 

which, on being substituted in the perturbative expansion of q, together with the change 
r — > t, gives us 



q (t) = ^-^^(o^-aiARATr/CW) + e ^t e 3iARAAr/(W) /3 t( )) . (4.5) 

We have thus obtained an asymptotic expression for this operator. On computing 
(n — l\q(t)\n), we see that the energy difference between levels comes out as E n — E n _\ = 
hu(l+ (3Xhn) I (4mw 3 ) + 0(A 2 )), consistent with all previous computations of this quantity. 

It has to be observed that our result is identical to the one obtained by Bender and 
Bettencourt 0, as is only to be expected, given the equivalence of the multiple scales method 
and the RG methods for a wide class of differential equations, to which the (classical) Duffing 
equation belongs. On the other hand, note the simplicity of our approach, where no a priori 
scale has to be assumed. 

In order to stress this latter point, let us consider the second order computation for this 
problem. The source term for q 2 is given by — {q$qi + <7o9i<Zo + <7i9o)j w here we have to 
consider the full q\ and not just the singular part. In this source term there will be terms 
that will give rise to secularities of the form e ±3lujt (t — r) and e ±tuJt (t — r) 2 . These we shall 
be able to ignore, because the renormalization to first order takes care of them. As a matter 
of fact, this is precisely what the statement of perturbative renormalizability amounts to in 
our case: that no divergences of a different form arise in the process of renormalization, that 
is to say, that all divergent (secular) terms can be taken care of by renormalization of the 
terms f3e~ iu)t and fte iwt . 

Another (simple) technicality in the problem at hand is that, since we have checked that 
is constant to order A 2 , we can use the commutator in the A and A 2 terms, thus 
making the computation somewhat easier. 

This results in the following expression for the secular relevant part of q 2 , q2,sr'- 

g 2 , sr = ~ 3 '!?7 r) ( (5A/" 2 - l)0t e ** - /3(5Af 2 - l) e -**) , (4.6) 



16a; 3 

whence the improved RG equation reads 



thus giving us 

and, as a consequence, 

E n - £ n _i = hu(l + (3A/m)/(4ma/) - 3AV(5n 2 - l)/(64mV) + 0(A 3 )) . (4.9) 



7 



V. QUANTUM PARAMETRIC RESONANCE 

As a last example of the usefulness of the RG method for quantum mechanical problems, 
we shall now illustrate its application to the phenomenon of quantum parametric resonance. 
Consider then the following hamiltonian: 

H = -^-P 2 + \muul(A + 2q cos(uj t))X 2 , (5.1) 

where A , q and u are constants. The evolution of any given state is computed by acting 
on it with the evolution operator U(t,t ), which satisfies 

ih^(t,t ) = H(t)U(t,to), (5.2) 

with U(t ,t ) = I, the identity operator. 

Let us divide the hamiltonian into an unperturbed and a perturbation part: 

H = H + H l= (^-P 2 + ^m^X 2 ) + (^mu 2 ((A - 1) + 2q cos(cu t))X 2 ) . (5.3) 

The reason for this decomposition lies in our previous knowledge that resonance will defi- 
nitely set in if A is equal to 1/4, but this is not essential for the final results. 
The evolution operator can be written as 

tf(Mo) = e-*-*> H °/ R EMMo) , (5.4) 
in such a way that the interaction picture evolution operator obeys the following equation: 

dJJ 

ih-^-(t,t ) = HjMUjfato) , (5.5) 

and in our case 

E l = ]^muol((A ~\)+ 2gcos(cc>ot))X 2 

= ^((A - J) + 2gcos(^ t)) (e-^-^a + e** ('-*>)/ V) 2 . (5.6) 

The constant operators a and at are the annihilation and creation operators at time to. 

We now perform the usual perturbative expansion for the interaction picture evolution 
operator, restricting ourselves to the Born formula, Uj — 1 — ^ J t * ds Hj(s). However, this 
leads to secular terms, and in order to eliminate them, we shall rather use this approximation 
close to the time t — r, by using the initial condition Ui(r,t ) = a(r), such that 

Ui(t,t ) = a(r) — — / ds Hi(s)a(r) + higher order terms. (5.7) 
Retaining only the secular terms, we obtain 

Uj(t,t ) = a(r) - ^(t - r) ({A - ^)(aa^ + at a) + q{e^ to a 2 + e - <woto (a t ) 2 )) a(r) + r. t. . 

(5.8) 
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For the sake of simplicity, let us set to = 0, without any loss of generality. We know that Uj 
cannot depend on the choice of r, and we are thus led to the RG equation to first order 



da iujQ 



{{A - i)(aa f + (Ja) + q{a 2 + (a f ) 2 )) a(r) = , (5.9) 



which, on being solved, provides us with an improved expression for the interaction pic- 
ture evolution operator, Uj(t,0) = exp(—itH e{i /h), where H e s is the large time asymptotic 
effective hamiltonian read off directly from the RG equation: 

H eS = ^ ((A - \){acJ + Ja) + q(a 2 + (a*) 2 )) 

= 2 (±(A - \ - q)P 2 + \rm%{A - \ + q)X 2 ) . (5.10) 

This is the first important result of our computation: we have resummed the effect of the 
variable frequency into an effective large time constant hamiltonian. If it happens that, for 
small positive q, \ — q < A < ^ + q, this effective hamiltonian corresponds to an inverted 
oscillator, thus marking the principal instability band (to the order we have computed). 

It now behooves us to compute the creation of particles due to this instability. In 
order to do this, we shall first write down the integral kernel that corresponds to Uj in the 



position representation, K(x, t; x', £'), using standard results for quadratic hamiltonians |T2 



Let 7 = uJ \Jq 2 ~ (A - 1/4) 2 and (p = J(q - 1/4 + A)/(q + 1/4 - A). The integral kernel 
K(x,t;x',0) := (x\Uj(t,0)\x') is computed to be (asymptotically): 

^■^•°) = (iSrt)^^(sTO^ + ^ )< ^ (7t) - 2ra ')) ' (5 ' n) 

It is now feasible to compute the asymptotic value of (n\U(t, 0) 1 0) through simple tabulated 
integrals, and we can calculate the transition probability from the ground state to even 
states: 

p a)= (2/)! 2v? ( V (512) 

U (/!) 2 2 2 7 V + (i + ^) 2s inh 2 ( 7 t) W 2 + (l + V 2 ) 2 ^\lt)J ' { ' ] 



It is easy to check that unitarity is preserved. An analogous computation leads to the rate 
of particle production, 

A^(t)= (1 ^ 2)2 sinh 2 ( 7 t). (5.13) 

These results can be compared with the computations of Shtanov et al. ||, and coincide com- 
pletely for the specific case at hand. Shtanov et al. arrive to this result through Bogolyubov 
transformations (to identify the function giving particle creation) and Krylov-Bogolyubov 
averaging (to perform the asymptotic analysis). This coincidence comes as no surprise given 
the first order equivalence of Krylov-Bogolyubov averaging to two-timing (a particular in- 
stance of the multiple-scales method) for a wide class of differential equations ||, and the 
(again first-order) equivalence of the multiple-scales and RG methods for many instances of 
equations. 
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VI. CONCLUSIONS 



We have performed several operator computations in quantum mechanics using the RG 
method for global asymptotic analysis. This method has the serious advantage that unitar- 
ity is built in, and that computations are simpler and more direct than in other techniques 
for asymptotic analysis. These ideas should be useful in a wide realm of applications. In 
the special case of quantum parametric resonance we derive explicitly an (asymptotic) ef- 
fective hamiltonian, which is an inverted harmonic oscillator whenever the system is in the 
instability region: the instability associated with the parametric resonance is turned into 
the unboundedness of the interaction hamiltonian, thus demonstrating the basic equivalence 
(asymptotically) of such different systems. Even so, unitarity is preserved throughout our 
computation, and the asymptotic results we obtain are well-behaved with respect to this 
fundamental property of quantum mechanical evolution. We have performed an analysis of 
the first instability band only for the quantum Mathieu equation. However, it is possible 
within this method to examine the whole forbidden/allowed band structure of this model, 
following in the quantum context the study carried out for the classical case by Goldenfeld 
and collaborators. As a matter of fact, for quadratic hamiltonians the whole instability anal- 
ysis can be reduced to classical mechanics, i.e., to classical Mathieu (or similar) equations. 
What we emphasise as novel in our results is the interpretation of instabilities as being due 
to effective large time inverted oscillator hamiltonians. Furthermore, analogous analyses 
can be carried for non-quadratic hamiltonians, even time dependent, where the quantum 
mechanical character of the problem would show itself to its fullest. 
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FIGURES 
FIG. 1. Backreaction: <p(t). 
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